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Divergence of the axial current and fermion density in Gross-Neveu models 
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The divergence of the axial current is used to relate the spatial derivative of the fermion density 
to the bare fermion mass and scalar/pseudoscalar condensates in 1+1 dimensional Gross-Neveu 
models. This serves as a useful test of known results, to explain simple features of the continuous 
chiral model and to resolve a conflict concerning the assignment of baryon number to certain multi- 
fermion bound states. 
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In this work we consider the family of Gross-Neveu 
models [l| with Lagrangians 

C = i#V - mo# + y [ W) 2 + A(^i75^) 2 ] (1) 

(A = 0, 1) in the large N limit. Flavour indices are sup- 
pressed, i.e., V>V = J2k=i^k^k etc. These are well- 
studied self- interacting fermion models in 1+1 dimen- 
sions, see 0, [H, H| and references therein. For vanishing 
bare mass mo they exhibit discrete (A = 0) or contin- 
uous (A = 1) chiral symmetry. For simplicity we shall 
refer to the first model as Gross-Neveu (GN) model, the 
second one as Nambu-Jona-Lasinio (NJL2) model in two 
dimensions. From the Euler-Lagrange equations we de- 
rive the divergence of the (flavor neutral) vector current 
j 1 * = ip^y^ip and axial vector current j§ = ipl^ls^^ 

V = 0, 

d,jg = 2 (m - (1 - \)g 2 ^) #75^ (2) 

The first equation expresses conservation of fermion num- 
ber in all models. The 2nd equation reflects the contin- 
uous chiral symmetry of the massless NJL2 model (for 
mo = 0, A = 1) and exhibits the source of chiral symme- 
try violation in those cases where the right hand side is 
non-zero. At mo 7^ 0, A = 1, it reduces to the standard 
PCAC relation. In 1+1 dimensions, vector and axial vec- 
tor current are trivially related (75 = 7°7 1 ), 
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35=3 



(3) 



Taking an expectation value of Eqs. f2J) in a static con- 
figuration (vacuum, multi-fermion bound state, dense 
matter) and using the factorization characteristic for the 
large N limit together with Eqs. (O, we find 

dilf) = 0, (4) 
Si (7°) = 2(m -(l-A) ff 2 (V^))(^i75V')- (5) 

Likewise we can take the thermal average of Eqs. @ 
in the grand canonical ensemble at finite temperature 



and chemical potential. Relations (|4I5|) then remain valid 
provided that the expectation values are interpreted as 
thermal ones. 

Whereas Eq. 0$ seems to be trivially fulfilled in all 
cases considered below and apparently contains little dy- 
namics, Eq. ([5]) turns out to be quite powerful. It relates 
the spatial derivative of the fermion density to the scalar 
and pseudoscalar condensates and the bare mass. In the 
remainder of this paper we will confront this identity to 
what is known about GN models to illustrate its poten- 
tial use. We will see that it provides us with a sensitive 
and non-trivial test of these models. In some cases, it 
sheds new light on simple physical properties which were 
known before. In one case it reveals an inconsistency in 
the literature and guides us towards its resolution. Vari- 
ous special cases of Eqs. (|4I5[) have been discussed in the 
literature before 0, [E| ■ 

Throughout this paper we use units such that the phys- 
ical fermion mass in the vacuum is m = 1. 

1. DHN kink-antikink baryons 

In the massless GN model with discrete chiral sym- 
metry ip —>■ 75^, [mo = 0, A = in Eq. (JTJ], Dashen, 
Hasslacher and Neveu (DHN) have constructed multi- 
fermion bound states by inverse scattering theory [||. 
These states can be viewed as Hartree-Fock (HF) solu- 
tions belonging to the reflectionless scalar potential 



S(x) = 1 + y(tanh£_ - tanh£ + ) 



with 



yx ± — artanhy. 



(6) 



(7) 
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The Dirac spectrum with potential ^ is symmetric un- 
der E — ► —E (charge conjugation) and exhibits two 
discrete states inside the mass gap at E = ± y/\ — y 2 1 
whereas the positive and negative energy continua are de- 
pleted by one state each as compared to the vacuum. The 
negative energy continuum is fully occupied (N fermions 
per level). Denoting the occupation numbers of the dis- 
crete states by n±, the self-consistency condition 

(W) = ~S(x) (8) 

g 2 

relates the parameter y to n± as follows, 

y = sm6, 9= ^-(n+-n_+JV). (9) 
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The parameter y controls the shape of S(x) as well as 
the mass of the baryon, Mb = 2Ny/ir. Let us now check 
Eqs. ([I} and ([5|). Using the fermion single particle wave 
functions from Ref. |7j we find that the expectation value 
of the current vanishes, (j 1 ) = 0, as pointed out already 
by Feinberg Q . This is consistent with Eq. ((4]) . Next we 
evaluate the fermion density (subtracting the density in 
the vacuum corresponding to S(x) = 1), 
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cosh £ + cosh 



where Nf is the fermion number 



N f 



n+ + ri- 



N. 



(10) 



(11) 



The pseudoscalar condensate in the GN model is rarely 
looked at, since it enters neither the self-consistency con- 
dition nor the calculation of the baryon mass. If we eval- 
uate it nevertheless in view of relation (0, we find 



1 



1 



(12) 



4 \ cosh 2 £ + cosh 2 , 

The axial current equation ^ reduces to 

d^f) =25(^i75^) (13) 

in the case at hand. (This relation was found earlier in 
Refs. 0,3.) Inserting the bilinears from Eqs. (|6I10I12|) 
into Eq. (1 13[) then yields the non-trivial identity 



yds 



\ cosh 2 £ + cosh 2 £ _ 
2[l + y(tanh£_ -tanh^ 



(14) 



1 



1 



cosh £ + cosh £_ 



which can indeed be verified for £± from Eq. ([7]). 

Knowing the pseudoscalar density is of interest for yet 
another reason. Any HF solution of the GN model with 
(xpi^ij}) = is also a valid HF solution of the NJL 2 
model. According to Eq. fT2")) . this is the case for all DHN 
bound states with Nf = 0, i.e., "baryonium" states. As 
can be seen from Eq. (jlip these states have a total of N 
fermions distributed over the two discrete levels (there 
are N + 1 such states, including the vacuum) and iden- 
tically vanishing fermion density. In the massless NJL2 
model, the axial current is conserved and Eq. §5§ reduces 
to di(j°) = 0. This is indeed satisfied by the GN bary- 
onium states with vanishing (j )■ These particular solu- 
tions of the NJL 2 model have already been discussed by 
Shei @. 

2. CCGZ kink baryons 

If y —> 1, the spatial extension of the DHN baryon be- 
comes infinite and it decouples into two separate entities. 
These can be regarded as another type of HF solution, 
the kink (or antikink) with scalar potential 



S(x) 



± tanhx. 



(15) 



This topologically non-trivial object is known in the lit- 
erature as Callan-Coleman-Gross-Zee (CCGZ) kink [y]. 
It has a single discrete state at zero energy which gets 
half of its strength from the positive and half from the 
negative energy continua. The fermion and pseudoscalar 
condensate are found to be (=F for kink and antikink), 



(i°> = Nf 

2 cosh 2 x 



N 



I 



2 cosh x 



(16) 



so that Eqs. (|5I13[) are again satisfied. The fermion num- 
ber has the value 



N 

N f =n . 



(17) 



The kink is a celebrated example for fractional fermion 
number @ (which we recover for odd N) and explains 
unusual spin-charge relations of solitons in polymers like 
polyacetylene [lol ] . Once again the baryonium state (n = 
N/2, Nf = 0) has vanishing pseudoscalar condensate and 
fermion density, solves the NJL2 model and is consistent 
with axial current conservation. 
3. Shei baryons 

The massless NJL2 model (mo = 0, A = 1) has contin- 
uous chiral symmetry ip — > exp(ia75)V'- Here the mean 
field approach involves a scalar and a pseudoscalar po- 
tential with two self-consistency conditions, 



; P(x). (18) 



Using the inverse scattering technique of DHN, Shei has 
identified massive baryons of this model Q . They were 
later confirmed by Feinberg and Zee within a resolvent 
analysis of the Dirac operator [ll|. Since the HF for- 
mulation of these bound states is not yet available but 
will be needed below, we have carried it out and included 
a brief summary in the appendix. In our notation, the 
self-consistent potentials of the Shei baryons are 



S(x) = l-/(a;)cos0, P(x) = f(x)shx9, 

,, . 2cos6» 
f(x) = 



1 + exp(2£) ' 



(19) 



with £ = x cos 9 and a parameter 9 living in the in- 
terval [7r/2,37r/2]. There is a single discrete state at 
Eq = — sin 9 which crosses the full mass gap as 9 is var- 
ied. Filling all negative energy continuum states com- 
pletely (N fermions) and the discrete state partially (n 
fermions), selfconsistency demands 



N 



(20) 



whereas the mass of the bound state can be shown to be 

N nn 

M B = — sin— . (21) 

7T N 
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Now we compute the fermion density, using the single 
particle wave functions from the appendix. First of all 
we find that the positive energy continuum is depleted 
by 1 — z^, the negative energy continuum by v in favor 
of the discrete state, with v = n/N. More specifically, 
the fermion density from the negative energy continuum 
(after vacuum subtraction) is given by 



<;°> 



2 cosh 2 £ 



(22) 



We recall that all continuum states are filled with N 
fermions. The appearance of the factor n is a combi- 
nation of the filling factor N and the depletion of the 
negative energy continuum states by n/N. The discrete 
state yields the density 



<i°)l 



discr 



n COS 6 

2 cosh 2 £' 



(23) 



where now n reflects directly the partial occupation of 
this state with n fermions. The total fermion density is 
the sum of both terms and vanishes identically, indepen- 
dently of n. The induced fermion density and the valence 
fermion density cancel exactly, and Nf — 0. Notice that 
the CCGZ kink with half filled discrete level is contained 
as a special case 9 = tt in the bound states found by 
Shei. As discussed above, it also has vanishing fermion 
number. 

Total screening of fermion charge as a result of a can- 
cellation between "valence" and "sea" fermions is at first 
sight surprising. Going back to Eq. (O, we recognize 
that the underlying reason is axial current conservation 
in the massless NJL2 model. It translates into the very 
restrictive condition di(j°) = for time independent con- 
figurations. Therefore, if the charge density vanishes at 
\x\ — > 00, it has to be identically zero. Chiral symme- 
try forbids a local accumulation of fermions. This ef- 
fect is qualitatively reminiscent of the depletion of soli- 
ton char ge by a finite chemical potential observed in 
Refs. fl2lll3j]. but it is more dramatic in the present case. 

It may be worth pointing out that the Dirac-HF Hamil- 
tonian with potentials (fH))) , after a global chiral rota- 
tion (see Eq. (j38|) from the appendix), coincides with the 
Hamiltonian employed in the context of irrational frac- 
tional fermion number [3, EH, EH- As shown in these 
works, the depletion of continuum states depends only on 
the asymptotics of the scalar potential and the value of 
the constant pseudoscalar potential. In Ref. [IB] dealing 
with relativistic field theory, the potentials were exter- 
nal fields without any self-consistency condition relating 
the occupation of the discrete level to the shape and size 
of the potential. In Ref. [Hi on diatomic polymers, self- 
consistency was only imposed on the scalar potential; the 
pseudoscalar one was a given (charge conjugation symme- 
try violating) parameter. Hence in these works screening 
of charge as discussed above was not an issue. 

In the present case, unlike in the GN model, we are in 
conflict with the original assignment of baryon number by 



Shei. The baryon mass is symmetric around the midpoint 
6 = tt of the interval in which 9 lives, cf. Eq. (|66|) of 
the appendix. Shei interprets this symmetry as charge 
conjugation. Filling the valence level with up to N/2 
fermions is supposed to describe baryons, from N/2 to 
iV" antibaryons. We believe that the correct assignment 
is that all of the states found by Shei have Nf = 0; they 
are no baryons at all, but rather baryonium states. A 
similar objection has been raised in the literature before 
by Farhi et al. [llj (see their appendix C), although no 
attempt was made to understand the vanishing baryon 
number. The same result can also be obtained from the 
resolvent analysis of the Dirac operator, see the Appendix 
A of § and Sect. Ill of [l8|]. Axial current conservation 
is a simple and strong argument in favor of these earlier 
works and the present study. The lesson one learns here is 
that one cannot read off baryon number from the filling 
fraction of the valence level(s) in cases where the self- 
consistent potential induces fermion number. 
4- Chiral spiral 

Another kind of baryon known in the NJL2 model is 
the massless one due exclusively to induced fermion num- 
ber ("chiral spiral" with unit winding number [l9|). Here 
it is preferable to work in a finite box of length L. The 
relevant potentials 

S(x)=cos , P(x)=sin (24) 

lie on the chiral circle S 2 + P 2 — 1, i.e., the vacuum 
manifold for spatially constant condensates. Since 



7°5(x) + i^Pix) = e^W^Og-is^W^ 



(25) 



one can map the HF problem with potentials ([24]) onto 
the vacuum problem with S(x) — 1, P(x) — by a local 
chiral rotation with linear ^-dependence. This induces 
spectral flow and increases the number of negative en- 
er gy s tates by B, as can be shown either by point splitting 
5, 19] or an appropriate cutoff method [2(|. The baryon 
density is constant (equal to B/L) so that the baryon is 
spread out over the whole space. The baryon mass van- 
ishes in the limit L — > 00. These massless baryons - also 
known from the 't Hooft model [U, - fit nicely into 
the present discussion, illustrating the fact that continu- 
ous chiral symmetry enforces spatially constant fermion 
density. So far, only integer baryon number was consid- 
ered (5 = 1 corresponds to N extra fermions, so that 
B = Nf/N). Presumably one could also allow for values 
Nt < N by letting the potential make only a fraction of 
a whole turn. This possibility was not envisaged so far, 
since it is most likely not relevant for dense matter and 
the phase diagram of the model. 

In summary, it seems that the chiral spiral type 
baryons are the only multifcrmion bound states of the 
massless NJL2 model which carry fermion number. They 
are completely delocalized as a result of axial current con- 
servation. 
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5. Phase diagram of massless NJL2 model 

In Ref. @ the phase diagram of the massless NJL2 
model at finite temperature T and chemical potential 
/j, was derived. The HF solution again exhibits a po- 
tential of helical shape (chiral spiral), with radius de- 
pending only on temperature and pitch depending only 
on chemical potential. In the present context, the most 
important finding is the fact that the fermion density 
is x-independent everywhere in the (T, fi) plane. The 
breakdown of translational invariance as a result of the 
Peierls instability (gap formation at the Fermi surface) 
only shows up in the (periodic) scalar and pseudoscalar 
condensates. These findings are consistent with Eq. ([5]) 
which predicts d\ (j ) = for the grand canonical ex- 
pectation value of the density. The reason behind the 
surprisingly simple phase diagram (as compared, e.g., to 
the GN model with discrete chiral symmetry [J]) then 
becomes clearer. 

6. Bare fermion mass corrections 

If we switch on the bare mass tuq in the GN model, 
Eq. (|13p remains valid with the scalar potential now given 
by S — mo — g 2 {tpip) . The discussion of the kink-antikink 
baryons can be repeated literally, the only change being 
the relation between y and the occupation numbers n±. 
The topological kink does not exist anymore because the 
degeneracy between the two vacua in the massless GN 
model is lifted by the bare mass term. 

Turning on mo in the NJL2 model is more interesting. 
Axial current conservation ([5]) gets replaced by the PCAC 
relation 



d^f) =2mo(^i 7 5^) = - 



2A 7 

7T 
-,2 



P, 



(26) 



where we have used (tpij^ip) = —P/g and expressed 
the ratio of the two bare parameters mo and g 2 by the 
physical parameter 



1 



It is related to the "pion" mass via [2c 



7 



: arctan ■ 



1 



n 



(27) 



A- (28) 



In Ref. [24[ the baryon in the massive NJL2 model has 
been determined approximately by means of the deriva- 
tive expansion, leading to a chiral perturbation series in 
powers of m 2 . To lowest order, the fermion density and 
the pseudoscalar condensate are given by 

U°) = -d xX , P=-sin2 X (29) 

where 2\ — 4 arctan e ni7 ' x is the sine-Gordon kink. In- 
serting Eqs. (f2"9"| into Eq. (j2"B"|) and using the leading or- 
der approximation 7 m^./4 to Eq. (|28p . we recover in- 
deed the static sine-Gordon equation for the chiral phase 
(£ = m^x), 



It was already obtained in Ref. [5[ and systematically im- 
proved by means of the derivative expansion in Ref. [24j . 
Turning now to the higher order corrections, we integrate 
Eq. (|2"6"|) to obtain the fermion density, 



(f{ x ))=- 2 -El I d.r'Pu'). 



(31) 



Another integration over x yields the fermion number 
Nf = N and the "sum rule" 



27 r°° 

1 = — - / dxxP(x) 
n J-00 



(32) 



sin(2x). 



(30) 



where we have integrated by parts and used P{— x) — 
—P(x). This relation provides us with a new way of 
testing the results of the derivative expansion. Insert- 
ing P(x) from Ref. [24| and 7 from Eq. (|2"5j) . we find 
indeed analytically that the right hand side of Eq. (|3Tj) 
is 1+0 (m®), consistent with the order of the calcula- 
tion. We have also checked that the higher order results 
of Urlichs [25| reduce the error to 0(m* 4 ), an excellent 
test of a lengthy and complicated analytical calculation. 
Relation (|26|) should also hold at finite temperature and 
chemical potential and may be of some help in deter- 
mining the (yet largely unknown) phase diagram of the 
massive NJL2 model |26| . 

Let us summarize our results. The most important 
equation of this paper is Eq. It relates the spatial 
derivative of the fermion density to the bare mass and 
scalar/pseudoscalar condensates in static configurations 
of GN models. As the derivation in a few lines shows, 
it is almost a triviality. Nevertheless it has many im- 
plications some of which have been exploited before 0]. 
The most striking predictions arise in the case where ax- 
ial current is conserved, i.e., the massless NJL2 model. 
Here the combination of axial current conservation and 
the close relationship between vector and axial currents 
in 1+1 dimensions constrains the theory most severely. 
As a result the fermion density is spatially constant in 
any static configuration of the massless NJL2 model, in 
agreement with previous calculations of massless baryons 
and the phase diagram based on the chiral spiral. Our 
result is in conflict with the Shei baryon, but this can be 
resolved by realizing that fermion density actually van- 
ishes as a result of a cancellation between valence and 
sea fermions. In those cases where the right hand side 
of Eq. ([S]) does not vanish, the identity has proven quite 
powerful for testing known results. We applied it suc- 
cessfully to the kink and kink-antikink baryons in the 
GN model, as well as to finite bare mass corrections to 
baryons in the massive NJL2 model. In all of these cases 
we confirmed the published results. Perhaps most note- 
worthy is the "sum rule" (f3"2"| relating the first moment 
of the pseudoscalar density to baryon number. It enabled 
us to test rather involved analytic calculations based on 
the derivative expansion to high order in a non-trivial 
way. 
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Appendix: Hartree-Fock treatment of Shei 
bound state 

Starting point is the Dirac-HF equation with the po- 
tential from Refs. [SGH, 

(75 jd x + 7 °(1 - f{x) cos 9) + i 7 7(x) sin ffj i/j = 

(33) 

where 

2k 



/(*) 



-, £ = k x, k — cos9. (34) 



1 + exp(20 : 
The global chiral transformation 

^ ^ e ie ^/ 2 ^ (35) 

maps the pseudoscalar potential onto a constant, 

(„1* +7 W-/M> + *-..)*-**. (36) 



Using the representation 



7° = -0"i, 7 1 = m, 75 = -0-2, 



the Hamiltonian now reads 

— sin 6* d x — cos 9 + f(x) 



H = 



(37) 



(38) 



-d x — cos 9 + f(x) sin0 
and becomes diagonal in the Dirac indices upon squaring, 





with 



11 ~ 1 -dl + U_ 



[/± = / 2 ±/' + l-2/cos< 



(39) 



(40) 



The resulting Schrodinger-type equations for the upper 
and lower components ip± of ip are identical to those of 
the kink in the GN model, 



cosh £ 



df + — ^ ; ) ^+ = -K 2 i>+ 



with 



E 2 -l 



K = ■ 

There is one discrete state, 

1 



1/2 



(41) 



(42) 



whereas normalized positive and negative energy contin- 
uum states are given by (E = ±\fk 2 + 1) 



4>e{x) 



1 



^2E{E - E 



ik — ko tanh £ 
E~E 



(44) 



We note the following contributions of single particle 
states to the condensates, 



V'oi75'0o 

4>Eils4>E 







1 



2 cosh 2 £ 
ko tanh £ 
E 
k 2 



(45) 



1 



2E cosh 2 £ l# - E E 
Self-consistency of the scalar potential 

^ 2 dk 



Eo 
E 



cos 6 — f{x) 



-A/2 27T 



(46) 



can be verified with the help of the vacuum gap equation 
(Ng 2 /ir) In A — 1. The self-consistency condition for the 
pseudoscalar potential, assuming n fermions in the dis- 
crete state, reads 

r A / 2 dk / \ t \ 

sm9 = -ng 2 <p i l5 i Jo -Ng 2 / — $ (47) 

Z7T 



A/2 



It holds for the following relation between 9 and v 

n/N, 



2 /l + sin( 

v = arctan 

7T \ cos 9 



(48) 



If v runs from to 1, 9 runs from 37r/2 to 7r/2 and -Eo 
from 1 to -1, thus the discrete level crosses the full mass 
gap from top to bottom. The fermion density for discrete 
and continuum states is 



Wo = 



2 cosh 2 i ' 



W e ^e = 1 - 



/■ 2 

h 



2E(E - Eo) cosh 2 £ 



(49) 



After subtracting the 1 in ip E tfjE, this was used to evalu- 
ate the fermion densities and continuum depletion factors 
in the main text. 

The solution discussed so far is related to the kink of 
the GN model. We now briefly sketch the construction of 
the antikink needed for evaluating the mass of the bound 
state below. Here, we start from the ansatz 

S(x) = 1 -f(x)cos6, P{x)=f(x)sm9 (50) 

where now 

f{x) = 2cos0-f(x). (51) 
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The relationship between / and / becomes more trans- 
parent upon writing 

f(x) — cos#(l — tanh£), 

f(x) = cos0(l+tanh£)- (52) 

Everything goes through as before except that / has to 
be replaced by /. In Eqs. (|39I41[) . this leads to an inter- 
change of upper and lower components. We now find the 
discrete state, 



cosh ^ 



E = sin 6, 



and the continuum states (E = ±\/k 2 + 1), 



1 



2E(E - E Q ) 



-(E-E ) 
ik — kg tanh £ 



(53) 



(54) 



The relation between the kink and antikink fermion wave 
functions is a discrete 75-transformation combined with 
the reflection 9 — ► 2n — 9. Assuming h fermions in the 
discrete state (occupation fraction v = h/N), condition 
gets replaced by 



2 /l-sin( 

v = arctan — 

7T \ cos 9 



(55) 



Notice that v + v = 1 for a common chiral angle 9 for 
kink and antikink. 

We are now in a position to evaluate the mass of the 
bound state. Since the kink and the vacuum cannot be 
described with the same boundary conditions (there is 
a chiral twist in the kink), it is not possible to apply 
directly the method from the GN model [6j, Ijj . Therefore 
we simply glue a kink and an antikink together, compute 
the mass of this whole object with periodic boundary 
conditions (the antikink un-does the chiral twist) and 
divide by two. 

For this purpose we need the sum over single particle 
energies from the discrete states, the HF double counting 
correction and the integral over (negative energy) contin- 
uum single particle energies. The discrete states give 

AA/sldiscr = nE + nEo = N(l - 2v) sin0. (56) 

The double counting correction is (using the gap equa- 
tion) 



AM 



2N 

7T 



dx- 



1 



cosh £ 
cos 9 In A. 



(57) 



For the continuum states we proceed as follows: Negative 
energy continuum spinors for the kink evolve as follows 
from large negative to large positive x [E^ — y/k 2 + 1), 



1 



y/2E k (E k -siat 
1 

y/2E k (E k -amt 



ik — cos 9 
-E k + sin 9 

ik + cos 9 
-E k + sin 9 



(58) 



For the antikink, the corresponding relation is 



E k + sin ( 



yj2E k {E k + sin 6») \ ik ~ cos0 

1 f E k + sin6» 

y / 2E k (E k + sin 6») V ik + cos0 



(59) 



When glueing the two objects together, we therefore have 
to multiply the antikink spinor by the phase factor 



C 



ik + cos 9 
^/k 2 + cos 2 9 ' 



(60) 



This maps the antikink spinor at large negative x onto 
the kink spinor at large positive x. We then find that the 
kink-antikink system produces the total phase shift 



5(h) = — arctan 



2k cos 9 
k 2 — cos 2 ( 



(61) 



This is common to upper and lower components, i.e., 
there is no chiral twist anymore. As a consequence the 
evaluation of AMs| cont is now straightforward [1,0] , 



AU i; |, m , = V / !|<W^ 



= AM B |surf + AM B |bulk- 

Using the asymptotics 



(62) 



^_2cos0 (fc ^ ±Qo) (63) 



and the derivative 



2 coa 9 



dk k 2 + cos 2 9 ' 

we find 

2N 

AM s | surf = COS0, 

7T 

2A" 

AMb bulk = cos 6» In A- NCI - 2i/)sin0 

7T 



(64) 



(65) 



The first part of AMs|bulk cancels the double counting 
correction (|57|) , the second part cancels the discrete state 
contribution (|56[) so that only AA/^lsurf survives. Divid- 
ing by 2 we finally get 



N 



N ■ N ■ - taa\ 

— aunrv = — sm7rj/ (66) 



in agreement with Refs. [8J, [llj 
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